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Abstract 

We analyze the prospects of the Tevatron for finding a Higgs boson in the two photon 
decay mode. We conclude that the Standard Model (SM) Higgs boson will likely not 
be discovered in this mode. However, we motivate several theories beyond the SM, 
including the MSSM, that predict a Higgs boson with enhanced branching fractions 
into photons, and calculate the luminosity needed to discover a general Higgs boson 
at the Tevatron by a two-photon invariant mass peak at large transverse momentum. 
We find that a high luminosity Tevatron will play a significant role in discovering or 
constraining these theories. 



Introduction 



In this letter, we investigate the possibihty that the Higgs boson(s) h associated with elec- 
troweak symmetry breaking may be found in the /i — > 77 decay channel at the Fermilab 
Tevatron. Our intention is to augment the many important studies preceding and associated 
with the Runll workshop on Supersymmetry and Higgs physics at the Tevatron |l| . In these 
studies, Standard Model Higgs boson detectability has been studied vigorously in its most 
promising production and decay channels. The classic channel of pp — Wh — * Ivhb remains 
the most important channel in the search for the SM Higgs boson, yet other modes can 
contribute to the total signal significance and perhaps yield evidence for the Higgs boson if 
sufficient luminosity is attained. 

We wish to study in detail the Higgs boson decays to two photons for many reasons. First, 
in our estimation this decay mode has not received adequate attention in previous studies. 
The capabilities of Higgs boson discovery in this mode should be carefully documented in 
order to better understand the Tevatron's full potential for Higgs boson detection. Second, 
there are many interesting and motivated theories that predict an enhanced decay rate into 
the 77 channel, and simultaneous suppression of the h ^ bb channel. Therefore, in these 
cases, non-standard search strategies must be employed to either find this Higgs boson or 
rule out its existence in the kinematically accessible mass range. And finally, we feel that 
studies such as these contribute to a more knowledgeable discussion regarding the worth of 
a higher luminosity Tevatron {e.g., run HI). 

Since there is no renormalizable and gauge invariant operator in the Standard Model that 
leads to — > 77 decays, it must be induced by electroweak symmetry breaking effects. The 
decay proceeds mainly through loop diagrams containing bosons and the t quark. The 
boson loop is dominant. The branching fraction for this decay in the light Higgs boson 
mass range 100 ^ rrih ^ 150 GeV is never much larger than 10"'^ since the 77 partial width 
must compete with the larger partial widths associated with bb, t~^t~, cc, gg, and WW* 
decays. The branching fractions for the Standard Model Higgs boson have been reliably 
calculated in Ref. 0. The maximum of the Standard Model branching fraction is 0.22% 
and is reached at ruh = 125 GeV. For rrih > 125 GeV, the branching fraction falls somewhat 
rapidly due to the increased importance of WW* decays. 
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Beyond the Standard Model 



It is widely recognized that the Standard Model is an unsatisfactory explanation of elec- 
troweak symmetry breaking. In this section, we review several well-motivated alternatives 
to the Standard Model Higgs sector. The first of these is low-scale super symmetry, where the 
symmetry breaking tasks are shared by two fields: and H^- receives a vacuum expec- 
tation value and gives mass to the up-type quarks, while receives a vacuum expectation 
value and gives mass to the down-type quarks. Both {Hu) and {Hd) vacuum expectation 
values contribute to the and Z° masses. In general, the sharing of the electroweak 
symmetry breaking task between two or more fields will disrupt expectations of Higgs boson 
phenomenology based solely on the analysis of the Standard Model Higgs boson. It is im- 
portant to identify regions of parameter space where our naive expectations fail, and where 
a more expansive search strategy must be engaged to find evidence of a Higgs boson. 

The mass matrix for the CP-even neutral Higgs bosons of supersymmetry in the {i?^, 
interaction basis is 

^ ( m\si\)? (3 + m\ cos^ (3 - sin (3 cos (3 {m\ + ni\) \ ^ ( Add A.^d \ q-, 
\ ~smPcosP{m\ + ml) m\ cos"^ P + m% sin'^ p J \ A^d ^uu J ' 

where m\ represents the pseudo-scalar mass, whose value is set by supersymmetry breaking, 
and Aij are quantum corrections whose form can be extracted from Ref. 0]. 
In the limit rriA ^ fnz the mass eigenstates of the above mass matrix are 

/i°ght = cos(3H', + sm(3H'^ (2) 
/iLvy = -sin/5i/° + cos/3i/°. (3) 

One can immediately see that (/^ught) = and (/iheavy) = 0) it is also true that all 
interactions of /^nght are equivalent to the SM Higgs boson. It is instructive to rotate the 
Higgs mass matrix to the {/^-ught; ^heavy} basis: 

,2 _ / m|cos2 2/5 -m| sin 2/5 cos 2/5 \ / A'^^ A[^ \ 
- y -m| sin 2/3 cos 2/5 + m| sin^ 2/3 )^ [ A'^^ A[^ ) ' ^ ' 

where the A-^ can be expressed in terms of the more commonly given corrections Aij, 
A[^ = AddCos^/3 + 2A„dCos/3sin/3 A„„sin^/3 
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A'la = -^dd cos /5 sin /? + A^^ cos 2/3 + Auu cos /3 sin /? 
A22 = Add sin^ /3 - 2Aud cos /3 sin /3 + A„„ cos^ /5. 

"Higgs decoupling" in supersymmetry means that one Higgs boson stays light and couples 
just like the SM Higgs boson as supersymmetry breaking mass scales get very high. This 
property of the supersymmetric Higgs sector can be immediately understood as a complete 
SU (2) Higgs doublet becoming very heavy (/iheavy^ H"^), while another doublet stays light 
(/i°ghj, W^). In the expressions above, this is equivalent to noting that m\ occurs only 
in the Al'22-element of the /i^ght ~ ^Lavy ^lass matrix. 

In supersymmetry model building, the supersymmetry breaking scale is a free parameter 
and is cycled over a very large range. This gives the false impression that over the vast 
majority of the parameter space, rriA is sufficiently larger than mz to be in the "decoupling 
region" described in the previous two paragraphs, and the lightest Higgs boson is well ap- 
proximated by /iilghf However, a natural electroweak potential — meaning a potential that 
has no large cancellations to produce the Z boson mass — prefers supersymmetry breaking 
near the weak scale. If we take naturalness and fine tuning arguments seriously, we expect 
tua ~ fnz, which leads to potentially significant deviations of the light Higgs boson couplings 
to the SM particles. 

An interesting departure from SM Higgs phenomenology occurs when the light Higgs bo- 
son mass eigenstate of supersymmetry is the weak eigenstate [Q-f^l- This scenario, or close 
approximations to it, can naturally occur in theories with large tan/5 = {Hu) / (Hd) , which 
are motivated by supersymmetric 5*0(10) unification f^, and by minimal gauge-mediated su- 
persymmetry theories that solve the soft CP- violating phase problem . The eigenstate 
has no tree- level coupling to bb or t^t^, and the total width for this light Higgs boson is 
greatly reduced. Loop corrections can modify these arguments. For example, supersymmetry 
breaking can induce couplings such as X'^H*bb (and X'^H*bb) in addition to the usual XhHdbb. 
The most important of these corrections often comes from gluino-squark loops (which do 
not contribute to A;^). 

If a significant \[ coupling is induced, the condition for shutting off the bb coupling is to 
shift the Higgs rotation angle to tana = e/ tan/3 (or tana = — e/ tan/5 for the case when 
the heavier CP even Higgs boson is Standard Model-like) where e = Ami,/ {nib ~ Arrih) and 
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Amfc = X[{Hu). The rr branching is not zero now, but is modified by a factor of compared 
to the Standard Model (in the hmit that A!^ is small). In contrast to the suppressed down- 
type fermion couplings, the partial width to two photons is equal to that of the Standard 
Model since no down-type quarks or leptons contribute significantly to the loop diagrams 
in either case. For these reasons, the branching ratio for the two photon final state can be 
greatly enhanced. Furthermore, the production rates through gg and qq' Wh^ are 

the same as in the Standard Model, since neither of these rely on the down-type fermion 
couplings. 

Other interesting theories imply enhanced branching fractions to two photons. The 
bosonic Higgs ^^jj that gives all fl^, or rather nearly all, the mass to the vector bosons 
and has no couplings to fermions is a good example of a Higgs boson with enhanced branch- 
ing fractions to two photons. However, the production cross-section of gg /i^^ is negligible 
in this model since the top quark does not couple to this Higgs. One must rely completely on 
electroweak boson couplings for production of the h^^i^ such as in qq' — >• Wh^^^i or WW h^^. 

Another example that has suppressed couplings to the fermions is an electroweak Higgs 
boson /ig^ added to top-quark condensate models [|ri|]. In this approach, the top and bottom 
quarks are assumed to get their masses through a strongly coupled group that condenses top 
quark pairs |T^, and all the remaining fermions and vector bosons get mass mainly through 
{/ig„). A good approximation in studying the phenomenology of a light /i^^ is to assume that 
it couples like the Standard Model Higgs to all particles except the top quark and bottom 
quark, to which it has zero couplings. 

In Fig.|l| we plot the branching fraction into two photons for the four Higgs bosons that 
we mentioned above: h^^, /i^j^, and h'^^. In each non-SM case considered, the branching 
fraction is larger than that of h'^^. Some models that we have not discussed here may have 
even higher branching fraction or perhaps lower. It should be kept in mind that any model 
of physics beyond the simple Standard Model will likely have different branching fractions 
into two photons. Since the two photon partial width is a one-loop process, it will also be 
sensitive to new particles in loop diagrams. Hence, even greater variability is possible than 
what we have shown here. For example, supersymmetric partners in the loops may increase 



or decrease the overall partial width of 77 |]T3|, |T4|, |T5| . In general, we should be prepared 
to discover and study a Higgs boson with any branching fraction to two photons, since that 
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Figure 1: Branching fraction into two photons for four different types of Higgs bosons: (1) 
h^^ is the Standard Model Higgs boson, (2) is the Higgs boson with Yukawa couphngs 
only with up-type fermions, which can be a mass eigenstate in large tan (5 supersymmetric 
theories, (3) is the Higgs boson that may help complete top quark condensation models 
as described in the text, and (4) h^-^ is a Higgs boson with tree level couplings only to W 
and Z bosons. 
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is perhaps the most hkely branching fraction to be altered significantly by new physics. 

There are several sizable sources of Higgs boson production within the Standard Model. 
At the Tevatron, they are gg — > /i, which is the largest, followed by qq' — Wh and qq 
Zh. For a heavy enough Higgs boson, the vector boson fusion processes WW, ZZ h 

2 

are also competitive. Although the decay branching fraction of — > 77 is of order yf^) 
there is hope that the narrow M^^ peak of the signal can be utilized to cut extraneous two 
photon backgrounds to sufficiently low levels that a signal can be detected. In the following 
two sections, we discuss search strategies based on inclusive and exclusive final states. A 
description of our calculational methods is provided in the Appendix. 



Inclusive 77 + X production 

First, we consider the total inclusive production of Higgs bosons, followed by their prompt 
decay to 77, where all Higgs boson production mechanisms can contribute. To study inclusive 
production, we apply cuts only on the properties of the individual photons or the photon 
pair, without studying the rest of the event in great detail. Before applying any cuts, the 



photon energy E'^ is smeared by a resolution function typical of Run I conditions pi 

AE^ .15 

- .03. (5) 



Je^ ( GeV) 



To optimize the acceptance of signal events, while reducing the "irreducible" backgrounds 
and those from jets fragmenting to photons, we apply the cuts, 

pTp > 20 GeV, \ri'~'\ < 2 (triggering and acceptance) 



Ai?^^ = ^{r]J - rilY + (07 - (fiY > 0.7 (separation) 

-p^<2GeV (isolation). (6) 

{i),R<A 

The high pt, central photons constitute a suitable trigger. We failed to find a more efficient 
Pj, cut than the one listed. With these cuts, the dominant source of background comes from 
the qq 77 process. For each event, we treat the diphoton pair as a Higgs boson candidate 
with mass M^^. 
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A further cut on the angular distribution of the photons in the rest frame of the Higgs 
boson candidate increases S/B: 



|cosr|<0.7. (7) 

The angle 9* is defined to be the angle that the photon makes with the boost direction in 
the 77 rest frame. The signal is rather flat in cos9* whereas the raw background peaks at 
I cos^^*| = 1. This cut is somewhat redundant to the other acceptance cuts, but will suppress 
fake backgrounds. 

With the above cuts it would require well over 100 fb^^ of integrated luminosity to even 
rule out a SM Higgs boson (95% C.L.) at any mass (the details will be given later). Therefore, 
new physics that provides a significant enhancement of the 77 + X total rate is required for 
this kind of signal process to be a relevant search. Large enhancements can occur either in 
the production cross-sections or in the decay branching fraction to photons. In the Standard 
Model, the gg h process constitutes roughly 2/3rds of the total production rate up to 
about m/j = 160 GeV, while the rest of the rate is mainly qq —>■ W/Z + h production. 
One does not expect the production cross-sections qq W/Z + h to ever greatly exceed 
the SM production cross-section given the nature of Higgs boson couplings to electroweak 
vector bosons. One does expect, however, that the gg h rate could be greatly enhanced 
by an increased coupling of the top quarks to Higgs boson |TB[, or by many virtual states 



contributing to the one-loop, effective ggh coupling, or from higher dimensional operators 
induced in theories with large extra dimensions |l^. In fact, these effects that increase the 
rate of gg h production will usually also alter the h ^ '-ff branching fraction. Therefore, 
we focus on the possibility of large enhancements of the ratio 

^ ^ (^(99 h)B{h -» 77) 

a{gg h)s,nB{h 77)sm ' 

In the following we will investigate the inclusive 77 + X rate from gg ^ h signal produc- 
tion alone. Although this underestimates the total cross-section by not taking into account 
the W/ Z + h and WW, ZZ — > h contributions, it lends itself to easy generalizations of 
Rgg ^ 1 where there is hope to find a signal at reasonable luminosity and where the other 
contributions are very small in comparison. Later on, we will see that the qq — W/Z + h 
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mh[GeV] a,ig(77 + X) [fb] da^^^J d M^^ [ih/GeV] 



100 1.49 39.3 

110 1.43 27.6 

120 1.29 20.2 

130 1.02 15.1 

140 0.73 12.3 

150 0.42 9.8 

160 0.13 7.4 

170 0.029 5.7 



Table 1: The total 77 + X production rate (in fb) for a Standard Model Higgs boson, and 
the differential rate (in fb per GeV) for the Standard Model backgrounds after applying cuts 
Eqs. (§-(0). 
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1.52 
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2.08 
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Table 2: The invariant mass resolution for a narrow signal from our simulations. The 
resolution is the la value of a Gaussian fit to the simulated signal after applying cuts Eqs. (|^)- 

m- 

signature alone lends itself to a useful, complementary analysis based on exclusive final 
states. In Table |l| we list the total SM signal cross-section {gg h only) and the differential 
background rate after all cuts have been applied. 

The Higgs boson width in the Standard Model is less than 20MeV for rrih < 150 GeV. 
Therefore, the invariant mass measurement of the two photons will have a spread entirely 
due to the photon energy resolution of the detector, which we call AM^^. In Table Q we 
show AM^^ for various Higgs boson masses, based on folding the photon energy resolution 
function Eq. (|^) with the photon kinematics. 

Based on Tables |1] and |^, we are now able to determine the significance of the signal with 
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respect to background after all cuts. We use the formula, 

S 0.96(TsigV^ 



where 



^^ = 71= • ^ (9) 

abkgd = 4AM,,^, (10) 

and C is the integrated luminosity. This formula counts the significance of signal to back- 
ground within a mass window M^^±2AM^^. This is a conservative and simple choice. When 
both the signal and background can be described adequately using Gaussian statistics, and 
the signal itself has a Gaussian shape, and the background is constant, the optimal mass 
window is M^^ ± v^AM^^. In our case, the background is not a constant, but the differ- 
ential distribution is well approximated by a straight line with a negative slope. Therefore, 
an asymmetric mass window (with respect to the peak) would most likely yield the best 
significance. We also require everywhere in our analysis that no limit or discovery capability 
is possible unless at least 5 events are present in this 2a spread mass bin. On the graphs 
we show below, this is a limitation mainly for the 2 fb~^ integrated luminosity curve. From 
Eqs. (P) and ([10|) , it is worth noting that an increase in integrated luminosity is equivalent 
to an improved energy resolution. 

In Fig. H we plot the 95% C.L. {Ns = 1.96) exclusion curves for a given luminosity in 
the Rgg-rrih plane. The SM Higgs boson corresponds to Rgg = 1 across the plot. We have 
put on a line on the graph corresponding to to give a non-SM reference example of Rgg. 
Other theories such as those discussed above can have Rgg much greater than that of h^^ or 
/iJJ. The plot is intended to be useful for comparing any theory to Tevatron capabilities. 

For a given integrated luminosity, the region above the corresponding curve can be ruled 
out to 95% confidence level. Therefore, with 30 fb^^ one could exclude a hi up to 120 GeV. 
The solid lines never cross Rgg = 1 which indicates that the SM Higgs boson could not be 
excluded in the 77 mode by the Tevatron even with over 100 fb^^ of data. One interesting 
limit to consider is a Higgs boson with only one-loop decays to gg and 77 final states. In this 
case, the production cross section x branching ratio is proportional to r{h gg)BR{h —>■ 
77) ~ r(/i ^ 77), and Rgg = r(fe,.,^oi)Bm,.M^^^) - 10'fSS&- Therefore, large values 



of Rgg are not unreasonable. 
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Figure 2: 95% C.L. luminosity contours of Higgs boson detection in the Rgg-rrih plane. For a 
given luminosity curve, the region below the curve cannot be ruled out with Tevatron data. 

Discovery of Higgs bosons with enhanced 77 + X production rates requires higher signifi- 
cance. For Ns = 5 we plot in Fig. ^the necessary enhancement Rgg to see a signal at this level 
at the Tevatron. With less than 30 fb~^, discovery is not likely for a h^^ Higgs boson with 
mass greater than 100 GeV. Therefore, Tevatron detection sensitivity in this channel is not 
as good as the Higgs boson search capacity at LEP2, which should exceed 105 GeV for both 
/igjjj and h^. Nevertheless other theories with larger enhancements of (j{gg h)B{h 77) 
may be discovered in the 77 + X mode first. 

Exclusive W/Z + h ^ W/Z + 77 signal of Higgs bosons 

We now attempt to gain more significance of signal to background by employing additional 
cuts. It is well known that the kinematics of resonance production at hadron colliders 
can be significantly affected by multiple soft gluon emission. Because of the different color 
factors associated with the qq 77 and gg ^ h processes, the pj^ spectrum of the Higgs 
boson signal is harder than the background. One strategy of LHC searches is to exploit this 
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Figure 3: 5a discovery contours of Higgs boson detection in the Rgg-ruh plane for a given 
luminosity. For each luminosity curve, the region above the curve can be discovered with 
Tevatron data. 



difference to establish a Higgs signal []TH|. However, the process W/Z + h ^ W/Z + 77 
typically has large pj^ even before including these QCD effects. At the Tevatron collider, 
the W/Z + h production process is relatively much more important than at the LHC, and 
quickly becomes the dominant process at even moderate values of pj^ with respect to M^^. 
For this reason, and because many of the extensions of the SM considered here have no 
gg ^ h component or only one of SM strength, we concentrate only on the W/Z + h signal 
in the following. For reasons discussed later, the WW/ZZ —>■ h signal is not as relevant for 
our analysis. 

We have done an analysis of varying the pj^ cut to maximize the total signal significance. 
We find that we optimally retain a significant portion of the total Higgs boson signal while 
reducing the backgrounds with the requirement, 

> (11) 

Also, the two photons from the Higgs boson decay tend to be balanced in pt, so we demand 

pI > (12) 
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as a further aid to reduce backgrounds and increase Sj B. 

After demanding such a significant cut on the dominant background becomes 77+I 
jet. However, the signal will most likely not have this topology. Rather, the decay of W 
and Z bosons can lead to: (1) two hard jets with Mjj ~ Myj/, M^, (2) one or more high pt 
leptons from W —>■ e, fi and Z — > ee, /i/x, or (3) missing transverse energy from Z uu, 
W T ^ soft jet, and W ^ e, fi ^ soft or very forward leptons. Therefore, it is useful 
to consider 77 signals that have one or two leptons, or missing energy, or two leading jets 
with Mjjj2 ~ mw,mz- To this end we require at least one of the following "vector boson 
acceptance" criteria to be satisfied: 

(a) > lOGeV and Ir/^'^^l < 2.0. 

(b) E/t > 20 GeV. 

(c) 2 or more jets with 50 GeV < Mj^j^ < 100 GeV. 

To perform this analysis, we constructed jets {Ej. > 15 GeV, \ri^ < 2.5 and R = 0.5) using 
the toy calorimeter simulation in PYTHIA with an energy resolution of 80%/ ^Je^{ GeV) . 
I/x was calculated by summing all calorimeter cells out to 77 = 4. 

The vector boson acceptance cuts eliminate a fair portion of the 77 plus jet background, 
as well as a potential contribution from vector boson fusion. The total rate of the vector 
fusion process (without cuts) is comparable to W/ Z + h only for > 160 GeV, where large 
values of B{h 77) are not well motivated (see, e.g., Fig. |l|). Nonetheless, we examined 
the effects of replacing cuts (a)-(c) by the requirement Mj^^j^ > 100 GeV to accept the jets 
associated with vector boson fusion: qq —>■ q'q'h. The results were not as promising as those 
based on cuts (a)-(c), and so we did not include a Mj-^j^ > 100 GeV acceptance cut in our 
analysis. 

In Table |] we show the signal and differential cross-section rates after all cuts, including 
the > and "vector boson acceptance" requirements. Table ^ can then be used to 

determine detectability of a Higgs boson given its mass and Ry- 

a{W/Z + h)B{h-^^^) 
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ruh [GeV] a,^{^^ + V) [fb] da^^^d/dM^^ [fb/GeV] 



100 0.134 0.102 

110 0.119 0.070 

120 0.097 0.047 

130 0.077 0.034 

140 0.049 0.026 

150 0.026 0.020 

160 0.0082 0.014 

170 0.0016 0.011 



Table 3: The total signal (in fb) for •yy + V, where V represents additional states passing the 
"vector boson acceptance" criteria enumerated in the text. The last column is the calculated 
background given the same cuts. 

The parameter Ry is useful if we make the reasonable assumption that increases in a{W + h) 
and a{Z + h) scale equivalently. 

In Fig. ^ we plot the 95% CL. {Ng = 1.96) exclusion curves for a given luminosity in 
the Rv-rrih plane. On the curve we have put lines for and the purely gauge coupled 
Higgs boson h^. The SM Higgs boson corresponds to Ry = 1 across the plot. For a given 
integrated luminosity, the region above the corresponding curve can be ruled out to 95% 
confidence level. The luminosity curves never cross Ry = 1 which indicates that the SM 
Higgs boson could not be excluded in the 77 mode by the Tevatron even with 100 fb"^ of 
data. However, with 30 fb~^ one could exclude h^-^ 137 GeV and h^^ up to 129 GeV in 

the 77 channel alone. 

For Ns = 5 discovery we plot in Fig. |^ the necessary enhancement of B{h —>■ 77) to see a 
signal at this level at the Tevatron. Discovery is possible up to 126 GeV for the bosonic Higgs 
boson as long as at least 30 fb~^ is obtained, and can be discovered up to approximately 
114 GeV. Both discovery reaches are beyond the expected reach capacity of LEPII. 
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Figure 4: 95% C.L. luminosity contours of Higgs boson detection in the Ry-mh plane. For a 
given luminosity curve, the region below the curve cannot be ruled out with Tevatron data. 
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Figure 5: 5a discovery contours of Higgs boson detection in the Ry-mh plane for a given 
luminosity. For each luminosity curve, the region above the curve can be discovered with 
Tevatron data. 
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Discussion and conclusion 



We have analyzed the capabihty of the Tevatron to find a Higgs boson decaying into two 
photons. We have found that the SM Higgs boson cannot be probed beyond LEP2 capabil- 
ities if the Tevatron accrues less than 100 fb~^. However, Higgs bosons in theories beyond 
the Standard Model may be probed (discovered or excluded) effectively with significantly 
less luminosity. For example, a Higgs boson that couples only to the vector bosons but 
has no couplings to the fermions can be probed up to 127 GeV with less than ~ 10 fb~^ 
integrated luminosity. In the MSSM, when h ~ h^, so that h ^ bb is suppressed and the 
W/Z + h{—>- bb) signal vanishes, our analysis shows coverage up to rrih = 114 GeV with 30 
fb~^, and exclusion capability up to rrih = 129 GeV. 

In an attempt to be as model independent as possible, we have presented graphs (Figs. ^ 
^) of exclusion and detectability as integrated luminosity contours in the plane of Higgs mass 
and Ri {Rgg and -Ry), where Ri parameterizes the enhancement of the 77 signal cross-section 
over the Standard Model. Therefore, any theories beyond the SM that have predictions for 
production cross-sections and decay widths of Higgs bosons can be compared with these 
graphs to attain an estimate of the Tevatron's capability. Of course, these figures are not 
applicable to a Higgs boson that has an intrinsic width greater than the detector resolution. 
The theories discussed above are far from this case. 

Finally, we comment on previous studies of 77 invariant mass signals at the Tevatron |ll9 |- 
25| . Much of the earlier work emphasized detectability at lower luminosities of h^T^ with 



m/jO^ < 100 GeV, where the branching fraction to two photons was 0{1). For example, this 
was the Higgs boson and the mass region covered in Ref. |jl9|. They also performed their 
simulations at the parton level and applied much looser cuts than our analysis, and much 



looser than those typically used in the experimental analyses of D0 and CDF in Run I |2T 



Another important analysis was completed very recently in Ref. [^] with results similar 
to ours, although the analysis differs in several ways {e.g., no cut). This study suggests 
that a signal of two photons combined with a single jet or two jets is an effective method 
to search for a Higgs boson in the two photon decay mode, but it did not utilize cuts as 
stringent as ours. A careful comparison needs to be made among all the observables in the 
various studies under precisely the same assumptions to ascertain which observables are the 
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most effective. And, of course, a combination of all useful observables should be employed 
to maximize our sensitivity to Higgs bosons. The suggested two photon observables outlined 
in this paper appear to be useful additions to the list of Higgs boson search observables. 

Acknowledgements: J.W. thanks Lawrence Berkeley National Laboratory for partial support 
in the participating guest program. SM thanks G.L. Kane for useful conversations. 



Appendix: estimation of signal and backgrounds 
Signal 

As mentioned in the text, we are mainly concerned with the processes gg ^ h and qq Vh 
(with V = W or Z) . The gg h ^ 'yy process is calculated based on 6-space resummation 
(see, e.g., Ref. [§6]), and performed to NLO accuracy. The total event rates (without cuts) 
agree with other fixed order calculations |^^. Since the multiple, soft gluon emissions are 
integrated out, the effect of isolation cuts must be determined by some other means. We use 
a constant isolation efficiency per photon ejso = 0.95 for these inclusive studies. Our results 
can be easily scaled if necessary to account for a different efficiency. 

The qq Vh process is calculated using PYTHIA [^, but multiplied by a constant 
i^-factor based on the resummed calculation of Ref. |^ . For completeness, the contribution 



of vector boson fusion processes were also calculated using PYTHIA without any effective- 
W approximation and no K factor. This process was never relevant for our analysis, for 
reasons discussed in the main text. 



Background 

The background estimate of the inclusive production of 77 pairs - where kinematic cuts are 
applied only on the properties of the individual photons or the diphoton pair - uses the next- 
to~leading order, resummed calculation of Ref. ||30l apphed to the 2.0 TeV collider energy. 
Since the resummed calculation integrates out the history of the soft gluon emission, the 
photon isolation efficiency must be estimated by another means, such as a showering Monte 
Carlo or from Z boson data. We use a constant isolation efficiency per photon ej^o = 0.95 as 
for the signal. No backgrounds from fragmentation photons (e.g. 7r°, t] — > 77) are included in 
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our numbers. The results of Ref. show good agreement with Run I data, and contain only 
a small component of fragmentation photons. For simplicity, we have ignored it entirely. Of 
course, the actual contribution from fragmentation photons depends critically on the isolation 
criteria and on the minimum pTp. Note that the resummed calculation does include the final 
state bremsstrahlung processes, e.g. qg — > 577. We find that our calculational method yields 



good agreement with the invariant mass distribution presented in Ref. without a large 
fragmentation component. 

To estimate the backgrounds toW or Z + 77, where the gauge bosons decay leptonically 
or hadronically, we need to determine the properties of the individual quarks and gluons 
emitted in the standard 77 production processes. This is not straightforward, since parton 
showering is accurate at describing event shapes but not event rates, whereas the hard NLO 
corrections to the 77 production rate are known to be important. For moderate values of pj^ 
relative to M^^, a fixed order (in as) calculation is as accurate in describing the kinematics 
of the photon pair as a resummed one (the transition between the two perturbative schemes 
is handled naturally in the resummation formalism, but the gluon emissions are integrated 
out). Therefore, we use the partonic subprocesses qq 77(7, qg 11 q, etc., to set the 
event rate, plus the parton showering method to reconstruct the full history of possibly 
multiple gluon emissions. For the gg 77+jets background, we use parton showering with 
the gg —>■ 77 process, but using the improvements of Ref. E^ to approximate the NLO 



corrections (the effect of using the exact pentagon diagram for the gg ng process is not 
important [^). The hard scale is set to the photon pair invariant mass. In all cases, we 
calculate the isolation efficiency explicitly. 
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